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Recent demonstrations of high-harmonic generation (HHG) at very high rep- 
etition frequencies (^100 MHz) may allow for the revolutionary transfer 
of frequency combs to the vacuum ultraviolet (VUV). This advance ne- 
cessitates unifying optical frequency comb technology with strong-field atomic 
physics 

aaaa. 

While strong-field studies of HHG have often focused on 
above-threshold harmonic generation (photon energy above the ionization po- 
tential), for VUV frequency combs an understanding of below-threshold har- 
monic orders and their generation process is crucial. Here we present a new and 
quantitative study of the harmonics 7-13 generated below and near the ioniza- 
tion threshold in xenon gas. We show multiple generation pathways for these 
harmonics that are manifested as on-axis interference in the harmonic yield. 
This discovery provides a new understanding of the strong-field, below-threshold 
dynamics under the influence of an atomic potential and allows us to quanti- 
tatively assess the achievable coherence of a VUV frequency comb generated 
through below threshold harmonics. We find that under reasonable experimen- 
tal conditions temporal coherence is maintained. As evidence we present the 
first explicit VUV frequency comb structure beyond the 3rd harmonic fli I2. fioj]. 

High repetition frequency HHG has recently been enabled through the use of femtosec- 
ond enhancement cavities jj} ^j. In addition to potentially transferring frequency comb 
techniques to the VUV, the enhancement cavity technique allows for greater harmonic pho- 
ton flux and a near-perfect Gaussian fundamental beam which provides for exceptionally 
clean and high signal-to-noise studies of the HHG process itself. In this Report, we utilize 
these techniques and present a new set of experimental and theoretical studies of harmonics 
7 through 13, generated below and near the ionization threshold in xenon gas with an in- 
tense 1070 nm laser pulse. These harmonic orders are of great interest for the development 
of VUV frequency combs and it is crucial to understand how fluctuations in the driving 
laser intensity influence the pulse-to-pulse coherence properties of the train of VUV pulses 
through the non-linear HHG process. 

Surprisingly, we find that harmonics as low as the 7th contain multiple contributions 
with different intensity-dependent phases. Our laser system permits a clean observation of 
interference between these phase contributions at unprecedented signal-to-noise ratios. In 
particular, we see experimentally a strong contribution from a component that has a large, 



2 



intensity-dependent phase and we show theoretically that it originates in semi-classical laser- 
driven continuum dynamics. We demonstrate that this phase is significantly larger than for 
above threshold harmonics due to the influence of the atomic potential on these low order 
harmonics. These measurements also provide a quantitative basis to evaluate the long-term 
temporal coherence of the VUV pulse train, a necessary condition for the establishment of a 
VUV frequency comb. Under reasonable experimental conditions, we find that the frequency 
comb structure should indeed emerge. To verify this, we demonstrate for the first time 
coherence between successive pulses in our 7th harmonic pulse train (A « 153 nm), which 

fin 

shows a coherence time ~5 orders of magnitude larger than that shown previously [111. 1121]. 

For this study, we utilize a passive optical cavity to enhance a high-power mode-locked 
femtosecond fiber laser at 1070 nm. We are able to achieve the peak power necessary for 
HHG at high repetition frequencies (136 MHz), which gives a very high signal to noise ratio 
and enables the pulse-to-pulse coherence measurements shown later. Figure 1 shows the 
experimentally measured on-axis yields of harmonics 7 through 13 as a function of the laser 
intensity in the center of the gas jet. Interestingly, all four harmonics exhibit complicated 
intensity-dependent yields, with steep increases interrupted by steps. For harmonics 7 and 
11 the intensity steps are more pronounced and occur at lower intensities than for harmonics 
9 and 13. The insets show experimentally measured far-field spatial profiles of each harmonic 
at an intensity of 2 x 10 13 W/cm 2 . At this intensity, harmonics 7 and 11 also exhibit strong 
off- axis halos in their spatial distributions whereas 9 and 13 do not. We also note that, in 
contrast to on-axis, the off-axis intensity-dependent yield displays no interference effects. 

The experimental results are very well reproduced by the theoretical calculations (shown 
in Fig. 1 for comparison) obtained via the coupled, non-adiabatic solutions of the time- 
dependent Schrodinger equation (TDSE) and the wave equation for a gas of xenon atoms 
exposed to an intense, 1070 nm laser pulse. Our approach is described in detail in [3], with 
the important difference that for the work described here we are directly integrating the 
TDSE numerically within the single active electron approximation [jjj]. We are thus treating 
the laser electric field and the atomic potential on an equal footing which is necessary to 
describe harmonics with photon energies below and close to the ionization threshold. As 

n 

initial conditions for the calculation we use the same parameters as the experiment [151 ]. 
The theory results reproduce the overall increase of the yield with intensity and exhibit 
prominent intensity-dependent steps in harmonics 7 and 11, and less pronounced steps in 
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harmonics 9 and 13. In addition, the positions of the steps are remarkably well reproduced 
by theory. 

Our experimental and theoretical results suggest that the observed intensity dependence 
is due to interference between contributions to the harmonic generation process which have 
different intensity-dependent phases. To explore these different contributions, we analyze 
the single atom intensity-dependent dipole moment d M) for each harmomc q in terms of its 
conjugate phase variable a, as described in detail in [lfjj]. From the calculated dipole phase 
we deduce the weight of each contribution, called a quantum path in the literature, over a 
range of intensities around Iq : via the transform: 



where W{I — 1$) is a window function centered on 7 [161 ]. The contribution characterized 
by phase-coefficient aj has a phase which is proportional to the intensity through (j)j = 
C£jU p (I) /huj, where U p = I /Auj 2 is the ponderomotive energy in atomic units, and uj is the 
laser frequency. For above-threshold harmonics, the phase coefficients obtained in this way 



correspond well to those predicted by the semi-classical model including the familiar 
short and long trajectories with phase coefficients of a\ « 0.2n and a 2 ~ 2n respectively fis| . 
To our knowledge, this analysis has not been applied to below-threshold harmonics before. 

Figure 2 shows the results of the quantum path analysis for harmonics 7—13 in xenon. 
All the harmonics exhibit multiple quantum path contributions; the two dominant ones 
have phase coefficients a « and a 2 ~ 2.5n — 3n. To understand the origin of the a 2 
contribution to the below-threshold harmonics, we have studied electron trajectories in a 
generalized semi-classical model in which the atomic potential is present. Crucially, we have 
found that in the presence of the atomic potential, low-energy electrons can lose enough 
energy to have less than zero total energy at x = (the position of the ion), leading to the 
emission of below-threshold harmonics. In the tunnel ionization model, this can only happen 
for electrons that have followed the long trajectory because they return to the position 
of ionization, at a time when the force from the laser field is against their motion. We find 
that the long trajectories accumulate an intensity-dependent phase that is larger than in 
the no-potential case, with a phase coefficient close to 3n for the range of intensities we are 
interested in, rather than the 2tt predicted with no potential. This is in good agreement with 
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the result in Fig. 2. There are no negative-energy returns for short-trajectory electrons in the 
tunneling model because they return to X{ at a time when the laser field further accelerates 
them toward x — 0. This leads us to designate the a contribution as a multiphoton process 
with no intensity-dependent phase, and not a generalized short trajectory. We note that 
harmonics 3 and 5 show only an a contribution. 

The calculations in Fig. 2 confirm our interpretation of the experimental results in Fig. 
1 as manifestations of the multiple generation mechanisms. The strong halos in the far- 
field spatial profiles of harmonics 7 and 11 are consistent with a large a 2 contribution to 
these harmonics at intensities below 2 x 10 13 W/cm 2 . This is because the large intensity- 
dependent phase, coupled with the Gaussian spatial profile of the fundamental beam, pro- 



duces a strongly curved harmonic wavefront [19f] . Although multiple generation pathways 
are familiar from HHG high above the ionization threshold, it has never been observed (or 
predicted) so far below threshold. Previous works have noted that harmonics close to thresh- 

n 

old can have large intensity dependent phases [20] , and that harmonics well below threshold 
can exhibit a step in the intensity-dependent yield, which was interpreted as an atomic 



resonance effect 



211 ]. Our results explain and extend these observations to harmonics well 



below and around threshold without invoking atomic resonances. Since it turns out that 
the intensity dependence of the dipole phase below threshold is due primarily to continuum 
dynamics, we have now a general quantitative connection between the laser intensity and 
the phase for these harmonics. 

This discovery allows us to connect these measurements to several other strong field phe- 
nomena observed in recent years. These experiments, which measured ionization dynamics 
rather than photon emission, all invoke motion in the combined field of the ion and laser field 



to explain the behavior of electrons with energies below the ionization threshold 
Our observations, and the high dynamic range associated with them, open the possibility 
to study the electron dynamics invoked in these experiments directly via the HHG process. 

Since the HHG process is extremely nonlinear, it has been a significant concern that 
frequency comb coherence could not be maintained through high harmonic generation. Pre- 
vious experiments have measured the coherence maintained in HHG but only over the time 
scale of a single pulse [ll, 12]. However, to produce a frequency comb, it is a prerequisite that 



coherence is maintained between successive pulses in the pulse train as this sets the frequency 
comb linewidth to be less than the repetition frequency. This comb linewidth determines 
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the ultimate frequency resolution for experiments using the VUV frequency comb [5] . The a 
parameters discussed above are a direct measure of the intensity-to-phase noise conversion 
in the below-threshold generation process that is directly relevant to the linewidth of a VUV 
frequency comb mode. For instance, using a parameters extracted from the experiment, we 
can estimate the percentage of the power that would be removed from individual 1 kHz 
linewidth comb modes for a laser source with a 1% RMS intensity fluctuations with a char- 
acteristic 1/f distribution for each of the generation processes |25|. The results are presented 
in Table 1. For comparison, we include a and phase values for harmonics produced high 
above the ionization threshold, for instance in argon at an intensity of 2 x 10 14 W/cm 2 . 

Our estimates of the intensity-to-phase noise put firm limits on the intensity noise that can 
be tolerated to preserve the comb structure in HHG, and it is clear from Table 1 that a very 
high degree of pulse-to-pulse coherence should be possible in principle. To put this question 
to rest, we demonstrate pulse-to-pulse coherence of the 7th harmonic. For this study, the 
harmonic radiation is collimated and sent through a cross-correlation interferometer with a 
delay set so that a pulse in the harmonic pulse train will interfere with the subsequent pulse, 
allowing us to measure the degree of pulse-to-pulse coherence [2fj]. Since the 7th harmonic 
radiation is absorbed in air and most solid materials, extreme care was taken in the design 
and construction of this interferometer with all necessary optical alignments being conducted 
under vacuum. 

Figure 3 shows the interference signal as the path length of the interferometer was scanned 
over 600 nm in about 2 seconds. The pulse-to-pulse coherence is clearly visible in the figure 
as a cross correlation that corresponds to the 7th harmonic wavelength of ~150 nm. The 
contrast ratio shown in the figure is only « 6%, however if we take into account imperfect 
collimation and severe power imbalance in the interferometer, we estimate that the maximum 
achievable contrast ratio was 10-15%. To provide a zero background for a significantly 
enhanced measurement contrast, we also used a lock-in detection by dithering the length 
of the interferometer at 392 kHz as the length of the interferometer was scanned. The 
demodulated signal is also shown in Fig. 3. 

Since imperfect alignment in the interferometer and pulse overlap can easily decrease the 
contrast ratio further, our measurement provides only a lower bound on the coherence. By 
measuring the pulse-to-pulse coherence we place a lower bound on the coherence time (<~ 10 



ns) which is roughly 5 orders of magnitude longer than that measured before (ill . Il2l ] . The 
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corresponding frequency comb linewidth is narrower than 20 MHz. 

In summary, we have found experimental and theoretical evidence for multiple contribu- 
tions to the generation of below-threshold harmonics. In particular, one of these contribu- 
tions is dominated by laser-driven continuum dynamics, in analogy with the semi-classical 
model for generation of high order harmonics. Our measurements of the spatial and spec- 
tral profiles and of the intensity dependence of the interference pattern all confirm that 
the intensity-dependent phase is larger for these low-order harmonics than for high-order 
harmonics due to the influence of the atomic potential. We used our measured phase coeffi- 
cients to estimate the coherence one can expect in both low-order and high-order harmonic 
frequency combs, and demonstrated that indeed, pulse-to-pulse coherence was maintained 
for the below-threshold harmonic generation. 



Methods 



The enhancement cavity 

To excite the enhancement cavity, we use a high-powered, frequency-controlled 1070 nm fiber 
laser system that supplies 100 fs pulses at 136 MHz with a pulse energy of 75 nJ [27|, l28l |. 
When the pulse train is resonant with the cavity (obtained by actively controlling the two 
independent degrees of freedom of the frequency comb, the repetition frequency and the 
offset frequency), we achieve an enhancement of 260, which translates to an intracavity 
pulse energy of 19 /xJ. Two 10-cm radius of curvature mirrors within the cavity produce a 
focal spot area of 960 /xm 2 that yields a peak intensity of 4 x 10 13 W/cm 2 . To generate 
harmonic radiation, we inject xenon gas near the intracavity focus using a glass nozzle 
with a 100 /im aperture and a backing pressure of 425 Torr. An XUV diffraction grating 
is utilized as one element of the enhancement cavity so that a portion of the harmonic 
radiation generated diffracts out of the cavity and impinges on a fluorescent plate 29] . By 
imaging the fluorescence onto a CCD camera, we are able to make measurements of the 
power and far field profile for individual harmonics. The exact location of the experimental 
laser focus relative to the gas jet, which is only known to within 200 /im, is inferred by 
aligning the experimental and theoretical result at the lowest intensity step for harmonic 7. 
This results in all the intensity step positions being remarkably consistent between theory 
and experiment. 
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Theoretical calculations 

In the semi-classical calculations, we performed a further check into the origin of the ~ 
phase contribution. We attempted to recover short trajectories with small phase coefficients 
and negative return energies by initiating electron trajectories with non-zero velocity at 
x = 0, simulating their ionization as a multiphoton process. Under the limited range of 
initial conditions where we do find such trajectories, their phase coefficients are also larger 
than in the no-potential case with a\ « 0.57T. Fig. 2 shows that such a\ contributions are 
not present in the fully quantum mechanical results. This confirms our conclusion that the 
a contribution is a multiphoton process. 

To extract values for the phase coefficients a and a 2 shown in Table 1, we made careful 
measurements of the inner and outer parts of the spatial (or spectral) profile of harmonic 
7 at a well known intensity [30|]. If we assume that the harmonic has a Gaussian spatial 
(temporal) profile with an extra intensity dependent phase, and a focal spot (duration) 
which is one half of the IR spot (pulse), we find a 2 ~ 2.7n and a < 0.2n. Likewise, using 
the intensity separation between the two prominent steps in Fig. 1 as the "period" of the 



3l|. All three 



interference between the two contributions, we extract a value of a 2 — a G « 3tt 
of these measurements are thus in excellent agreement with the theoretical predictions. 
Coherence measurements 

For the pulse-to-pulse coherence measurements, the harmonic radiation outcoupled with 
the intracavity diffraction grating was sent through a VUV interferometer (see Fig. 4). A 
MgF 2 beamsplitter, calcium fluoride lenses and aluminum mirrors optimized for operation 
in the VUV were used for the construction of the interferometer. One end-mirror of the 
interferometer was mounted on a piezoelectric transducer (PZT) so that the optical delay 
of the interferometer could be scanned and quickly modulated. To detect radiation at the 
output port, we used a conventional photomultiplier tube (PMT) with a thin layer of sodium 
salicylate on the front surface. The sodium salicylate fluoresces at a wavelength of around 
420 nm when struck by the 7th harmonic radiation, which matches very well to the peak 
detection efficiency of the PMT. To observe the interference with a better contrast, we 
applied a 392 kHz modulation to the PZT. The PZT and mounted mirror had a natural 
resonance close to this frequency allowing enhanced mirror travel. The demodulated signal 
is basically the first-order derivative of the DC-based interference measurement. 
* Corresponding author: Ye@JILA.Colorado.edu 
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TABLE I: Frequency comb parameters for below threshold harmonics {Tiu < I p ) at an intensity of 
I = 2 x 10 13 W/cm 2 and for above threshold harmonics (hu > I p ) at / = 2 x 10 14 W/cm 2 . Also 
listed is the rms phase noise (0 rms ) and an estimate of the power left in a 1 kHz carrier for an 
individual frequency comb mode with 1% integrated intensity noise in frequencies above 1 kHz in 
a 1/f distribution. 
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FIG. 1: Harmonic yield plotted as a function of the intensity at the center of the xenon jet. (a)-(d) 
show harmonics 7 through 13, respectively. The red filled circles are measured and the blue solid 
lines are theory (see text), both for the on-axis yield. Theoretical harmonic yields were scaled 
in magnitude to offer a clear comparison. The insets show the far-field spatial profile for each 
harmonic. Halos on the beam profile and strong oscillations in the intensity dependent yield are 
clearly visible in harmonics 7 and 11, providing clear evidence of multiple generation pathways 
with distinct intensity dependent phases. 
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FIG. 2: Quantum path distributions calculated for harmonics 7-13. The color scale has been 
normalized for each intensity and shows only relative strengths. Harmonics 7 and 11 show large 
quantum path contributions at a « and a « 3tt that manifest themselves in the spatial profiles 
with halos and intensity-dependent yields with large oscillations, as shown in Fig. 1. Harmonics 
9 and 13 only show prominent contributions at a « for intensities below ~ 3 x 10 13 W/cm 2 , so 
the features described for harmonics 7 and 11 are largely absent. 
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FIG. 3: 7th harmonic pulse-to-pulse coherence measurement, a) Pulse-to-pulse measurement setup: 
VUV interferometer with the delay set for pulse n in the 7th harmonic pulse train to interfere with 
pulse n+1. b) Figure depicting the electric field of pulse n in the 7th harmonic pulse train that has 
been sent through an optical delay line so that it interferes with the n+1 pulse in the train, which 
is shown as the dotted red trace. A cross-correlation signal results as the length of the delay line 
is scanned, which is shown in c). d) The blue closed circles show our experimental measurement 
of the cross correlation as the path length of the interferometer is linearly scanned, showing a 
contrast ratio of « 6% limited by the spatial mode overlap and the power imbalance between the 
two interfering pulses. The red open circles display the signal achieved from a lock-in detection of 
a 392 kHz modulation of the delay line. Both sets of data are fit with sinusoidal functions. Point 1 
in c) and d) shows the positions of a maximum negative lock-in signal where the cross-correlation 
signal at DC has a maximum negative slope. Point 2 shows the position of a zero lock-in signal 
with the DC cross-correlation at its maximum. 
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FIG. 4: Experimental setup for the demonstration of VUV pulse-to-pulse coherence. A cavity- 
enhanced Yb fiber frequency comb is used to generate harmonics in xenon gas. Photodiodes 
(pd) are used to monitor the reflected and transmitted fundamental light and maintain resonance 
between the frequency comb and the enhancement cavity [l|, 0]. An intracavity VUV diffraction 
grating (vuv-dg) is used to outcouple the 7th harmonic radiation, which is then sent through an 
cross-correlation interferometer. One mirror in the delay arm of the interferometer is mounted on 
a piezoelectric transducer (pzt), allowing scans of the differential path length. The 7th harmonic 
radiation at the output of the interferometer is detected by monitoring the fluorescence from a 
sodium salicylate coated plate with a photomultiplier tube (pmt). 
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